Mode hopping strongly affects observability of dynamical instability in optical 

parametric oscillators 
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Theoretical investigations of dynamical behavior in optical parametric oscillators (OPO) have 
generally assumed that the cavity detunings of the interacting fields are controllable parameters. 
However, OPOs are known to experience mode hops, where the system jumps to the mode of 
lowest cavity detuning. We note that this phenomenon significantly limits the range of accessible 
detunings and thus may prevent instabilities predicted to occur above a minimum detuning from 
being evidenced experimentally. As a simple example among a number of instability mechanisms 
possibly affected by this limitation, we discuss the Hopf bifurcation leading to periodic behavior 
in the monomode mean-field model of a triply resonant OPO and show that it probably can be 
observed only in very specific setups. 

PACS numbers: 42.65.-k, 42.65.Sf, 42.65.Yj 



I. INTRODUCTION 

Continuous-wave optical parametric oscillators 
(OPOs) are tunable sources of coherent light that have 
proved extremely useful in quantum optics or high reso- 
lution spectroscopy 0. They also have attracted great 
interest as model systems in nonlinear dynamics because 
they are based on the simplest optical nonlincarity, 
three-wave mixing, and are expected to exhibit complex 
dynamical behavior in some regions of parameter space. 
Indeed, many theoretical studies have predicted a variety 
of complex temporal and spatio-temporal dynamics (see, 

e.g., HUBHIEHBH In p articmar > tne 

simplest model of a triply resonant OPO (TROPO), 
the degenerate monomode mean-field model, was shown 
twenty-five years ago to display a Hopf bifurcation 
leading to periodic behavior 0, y( , and later to exhibit 
deterministic chaos Quite surprisingly, this in- 

stability has to our knowledge not yet been observed 
experimentally. Although oscillatory behaviors have 
been reported in several experiments [Til Il3l Il4| . 
they have been shown to stem either from thermal 
effects E2> 113 or from the interaction of transverse 
modes |l4j . Very recently, chaotic behavior has been 
evidenced in a TROPO 151 but is believed to be also 
linked to multimode operation. 

Most of the instabilities that have been investigated 
theoretically require large values of the cavity detunings 
and of the pumping rate. As the latter is obviously lim- 
ited by the power of the pump laser, it would only be nat- 
ural to identify it as the limiting factor when searching 
for the instability [T^j. As f° r detunings, they are gener- 
ally considered as parameters that can be easily tuned so 
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as to pull the system away from equilibrium and observe 
complex behaviors. However, there are two limitations 
on the values cavity detunings can take. The first is 
simply that the OPO threshold increases quadratically 
with detuning, and thus that operation is restricted to 
detuning values for which threshold power remains be- 
low available pump power. The second arises because of 
a phenomenon known as mode hopping: OPOs sponta- 
neously choose their operation frequency so as to operate 
on the cavity mode with smallest detuning. As the OPO 
is pulled away from a cavity resonance in order to increase 
detuning, a more favorable operating point appears near 
another cavity resonance and the system jumps to it. Be- 
cause mode hops in OPOs occur for variations of cavity 
length as small as a few nanometers, they make it diffi- 
cult to stabilize doubly or tripl y re sonant OPO and to 
achieve smooth tuning 0, [PA Il8|. Here, our concern 
is that mode hops prevent the OPO from achieving high 
values of the detunings and thus restraint the parameter 
range that can be explored. Such a limitation might very 
well preclude the experimental observation of dynamical 
behaviors predicted theoretically. 

In this paper, we discuss this problem in a simple ex- 
ample, the Hopf bifurcation that leads to periodic behav- 
ior in the monomode mean-field model of a TROPO. We 
show that mode selection in OPOs indeed prevents exper- 
imental observation of this instability in most practicable 
configurations. We first recall the basic properties of the 
monomode mean-field TROPO model and the conditions 
under which the Hopf bifurcation occurs. In the limit of 
infinite pump power, a simple lower bound for signal de- 
tunings at which the Hopf bifurcation can occur is easily 
obtained, which confirms that this bifurcation requires 
high values of the signal detuning. In a second part we 
derive the expression of the maximal value that signal 
detuning can reach before a mode hop occurs, which de- 
pends on the length and finesse of the cavity as well as on 
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crystal birefringence. By comparing the two bounds and 
searching for parameter regions where they are compati- 
ble, we find that mode hops generally keep the TROPO 
away from parameter ranges where the Hopf bifurcation 
can be observed, unless very high-finesse and very short 
cavities are used, which would make operation extremely 
difficult. This is confirmed by a numerical exploration at 
finite pump power of the parameter space of this model 
for various values of the cavity finesses. It shows that 
mode hopping rather than pump power is the limiting 
factor in order to reach instability. This result provides 
a plausible explanation of the fact that the Hopf bifurca- 
tion of the monomode TROPO has not yet been observed 
experimentally. It also calls for further investigations in 
order to determine whether mode hopping also interferes 
with other predicted instabilities. 



II. TROPO MONOMODE MEAN-FIELD 
DEGENERATE MODEL: HOPF INSTABILITY 

We now recall the main features of the simplest 
TROPO model, the degenerate longitudinally and trans- 
versely monomode mean-field model 0, Q . Light gen- 
eration in an optical parametric oscillator is based on 
parametric down-conversion in a nonlinear crystal of a 
pump photon into two lower-frequency photons called 
signal and idler. In a TROPO, the optical cavity en- 
closing the crystal is resonant for all three fields so as to 
minimize operation threshold. In the mean-field (a.k.a. 
uniform-field) approximation, the time evolution of the 
normalized amplitudes A s , Ai and A p of the signal, idler 
and pump fields inside the cavitycan be described by the 



following differential equations 
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A s 
Ai 



-(l + iA s )A s 
-(X + iAi)Ai 



A* A 



A* A 



Ap = 7 [-(1 + iA p )Ap - A 3 Ai + E] 



(la) 
(lb) 
(1c) 



where A s , Aj and A p are the detunings between the 
optical frequency and the frequency of the closest cavity 
resonance for the signal, idler and pump fields, respec- 
tively, and E is the pumping rate. The time unit is the 
cavity decay time of the signal field and 7 is the cavity 
decay rate for the pump. In this paper we focus on the 
stationary regimes of Eqs and are interested in de- 
termining when they become unstable to give birth to 
periodic oscillations depending on values of control pa- 
rameters. 

A little known property of Eqs. Q is that although the 
signal and idler fields are in principle distinct, their time 
evolution can be described by a single amplitude after 
transients have died out. First it should be noted that 
stationary solutions of Eqs. exist only for A s = Ai, a 
relation that can be shown to follow from photon num- 
ber conservation |19| . Replicating a similar calculation 
carried out in the analysis of a bimode model |20(, it is 



then easy to show that Eqs. I|la(l and Ijlbfl imply that af- 
ter a sufficiently long time, the amplitudes A s and Ai are 
equal up to a constant phase difference that can always 
be removed by a redefinition of the amplitudes. This is 
obviously linked to the fact that signal and idler pho- 
tons are twin photons created in the same quantum pro- 
cess. Without loss of generality, the asymptotic dynam- 
ics of the TROPO can then be modeled by the following 
normalized equations (degenerate model) describing the 
time evolution of the complex amplitude of the signal 
field A s and of the pump field A p jj] : 



A s = -(l + iA s )A s + A*A 



pi 

1 [-{l + iA p )A p -Al + E] 



(2a) 
(2b) 



For pumping rates above the parametric emission thresh- 
old given by 



2& = (1 + A»)(1 + A2), 



(3) 



equations J2J have non-zero stationary solutions which 
have been shown to fit accurately experimental observa- 
tions in the vicinity of threshold [ll| . When pump rate is 
increased, these stationary solutions can become unsta- 
ble through a Hopf bifurcation giving rise to oscillatory 

behavforUHH. 

tion is 
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A necessary condition for this bifurca- 
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which ensures that there is a finite pump rate Eh > E t h 
at which the stationary nonzero solution bifurcates to a 
periodic solution, which is is given by Q 



E% = 



7 2 (l + A 2 ) + 4(l + 7 ) 
-2(1 + 7) 2 [1 
(A p + A s ) 2 , 



-1 2 



2(1+A P A 3 ) 1 
7(1+A£) J 



(A p A s 



(5) 



At higher pump rates, the limit cycle born in the Hopf 
bifurcation undergoes a period-doubling cascade leading 
to chaos 0). 

As inequality is closer and closer to equality, the 
Hopf threshold Efj given by 10 becomes larger and larger 
and is eventually rejected to infinity. For a given max- 
imal pump rate available, E maxi whether the instability 
can be observed at fixed detunings inside the parameter 
domain delimited by (£Q) is determined by the inequal- 
ity -Emax > Eh(A p , A s ). Since we are interested here in 
specifying the unstable region by simple bounds on the 
detunings, we first assume that we have infinite pump 
power available. Under this approximation, whose valid- 
ity will be checked in numerical simulations in Sec. IIVI 
the Hopf instability domain is solely determined by in- 
equality 

The instability domains in the (A p , A s ) plane have 
been plotted in Fig. for different values of 7, and are 
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FIG. 1: Map (A p , A s ) for different values of the parameter 
7. The curves delimits the areas for which condition gj is 
fulfilled. 



operating point becomes more favorable. As we discuss 
below, this process tends to minimize frequency detuning 
and consequently to limit its maximum value, which can 
be derived analytically. 

As with any optical oscillator where an amplifying 
medium is enclosed inside an optical cavity, there are 
two main constraints that determine the operating con- 
ditions of an OPO: the first one is that the gain in the 
amplifying medium must be sufficient to overcome cavity 
losses over each round trip in the cavity, the second one 
is that the generated field must be nearly resonant with 
one of the cavity modes so that amplification by the gain 
medium is cumulative over successive round trips. 

In optical parametric oscillators, the nonlinear inter- 
action is optimal when the relative phases of the three 
interacting waves remain fixed during propagation. The 
phase-matching condition is easily formulated in a cor- 
puscular point of view: the conversion of one pump pho- 
ton into signal and idler photons must satisfy energy and 
momentum conservation: 



seen to be bounded away from the origin. Along their 
boundaries, J3J is an equality and defines the signal de- 
tuning as a function of pump detuning. It is easily found 
that the minimum absolute value that the signal detun- 
ing can take on the boundary, and hence in the instability 
domain, is: 

A^ in =mm{|A s |} = Vt(tT2), (6) 

and is obtained for |A p | = y/(2 + 7) Ay. Note that A^ in 
is roughly proportional to 7 and thus increases with it. 

Since still holds at finite pump power, the lower 
bound provided by I© is always valid. This clearly shows 
that the Hopf bifurcation of the monomode degenerate 
mean-field model can only occur for sufficiently high sig- 
nal detunings. In the next Section, we describe the pro- 
cess of mode hopping and how it limits the values that 
signal detunings can take leading to a maximal absolute 
value for the detunings. Whether the two constraints 
can be satisfied simultaneously will eventually determine 
whether the Hopf bifurcation can be observed experimen- 
tally. 

III. LIMITATION OF SIGNAL DETUNING DUE 
TO MODE HOPPING 

A. Theory 

In the theoretical analysis of Sec [HJ it was assumed 
that the frequency detuning of the signal field is a fixed 
parameter. This is however not entirely true, because 
the operating frequency of an OPO is not actually cho- 
sen by the experimentalist but results from a complex 
mode selection process. While small variations of the 
cavity length most often modify frequency and hence de- 
tuning gradually, sudden jumps will occur when a remote 



uj p = lu s + uji (7a) 
k p = k s + ki (7b) 

where uj p ^ s ,i and k p . s ^ denote the frequencies and the 
wavevectors of the pump, signal and idler fields. 

The pump properties being fixed, relations gen- 
erally single out unique values for the signal and idler 
frequencies and much of OPO design consists in ensur- 
ing that these values fall in the desired frequency range. 
However, while energy conservation l|7al) strictly holds for 
continuous- wave OPOs, momentum conservation may be 
satisfied only approximately because of the finite size of 
the nonlinear crystal. Thus there is a small but finite 
frequency domain around the exact phase-matching fre- 
quency where the nonlinear gain is adequate for OPO 
operation. It is usually the case that many cavity reso- 
nances fall inside this domain, and an important problem 
is to determine the resonances near which oscillation will 
actually occur. 

Monomode optical parametric oscillators behave as ho- 
mogenously broadened lasers in that when there are sev- 
eral frequencies for which gain overcome losses in the 
empty cavity, the one with the lowest threshold takes 
over by saturating the gain in such a way that compet- 
ing modes remain below threshold. Mode selection thus 
amounts to determining which operating modes have low- 
est threshold. If variation of raw gain in the domain 
around perfect phase-matching is neglected, and consid- 
ering pump detuning as a fixed parameter, this is equiva- 
lent to finding the allowed operating frequencies for signal 
and idler fields that have lowest frequency detuning, as 
expression (J3J) for parametric threshold shows. 

In doubly and triply resonant OPOs, this problem is 
made difficult by the fact that the signal and idler fields 
must be simultaneously resonant. Thus operation can 
only take place at coincidences between two frequency 
combs. In general the two combs have different periods, 
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either because the signal and idler fields have very differ- 
ent frequencies or because they are polarized along dif- 
ferent axes of a birefringent crystal (the so-called type-II 
phase-matching pHl. A singular configuration that we 
will not consider here is when signal and idler have identi- 
cal polarizations and frequencies (type-I phase matching 
at degeneracy). As a result, frequency tuning in OPOs is 
a complicated problem that has been studied very care- 
fully both in the type-I and type-II cases [TgL Il7| . 

The detunings of the signal and the idler fields are 
given respectively by Alo s = u> s — u>f and Au>i = uii—uf', 
where lo s and ioi arc the optical frequencies of the signal 
and idler fields and uf and tof 1 are the closest resonance 
frequencies of the cavity. For the sake of simplicity, we 
assume in the following that lo s and u>i are close to a 1:1 
ratio but the argument can be generalized easily to any 
rational number. 

Taking energy conservation (|7a[) into account, the total 
frequency mismatch Alo = Alu s + Aw, for the signal-idler 
mode pair is given by : 



Aw = 



(8) 



which does not depend of the specific oscillation frequen- 
cies lo s and LUi chosen by the system, but only of the pump 
frequency and the resonance frequencies of the cavity for 
the signal and the idler modes. For our purposes, Auj is 
the relevant quantity to consider as the individual detun- 
ings Alu s and Awj are proportional to it in a stationary 
state, as a consequence of energy and photon number 
conservation [l^. The mode pair selected by the OPO 
will be the one that minimizes the total frequency mis- 
match Alu so as to minimize threshold. 

Taking into account that uif and tuf- belong to two 
frequency combs specified by the free spectral ranges Slo s 
and Sujf. 



N a Su, 



Ni6u>i 



(9) 



where N s and iVj are the mode indices, and having in 
mind that \5uj s — 5u>i\ <C Slj s , 5u>i, we rewrite J5J as 



Alu = lo p ~ N(Slo s + 6<Ji) - AN(Slo s - 5u>i) (10) 

with N = (N s + Ni)/2 and AN = (N s - N t )/2. Under 
our hypotheses, we have AN <C N. Hence the last term 
in (|10|) is in first approximation negligible compared to 
the second term, and the integer value of N that min- 
imizes Alo is determined independently of AN. Then 
the optimization problem can be refined by searching for 
the integer value of AN that minimizes (|1U|) at fixed N. 
Studies of tuning properties of double or triply resonant 
OPO have shown that the operating modes are grouped 
into clusters, each cluster consisting of a sequence of ad- 
jacent modes 0,^). In (|TU|l . indicates the cluster 
and AN distinguishes between adjacent modes inside the 
cluster. 

An important consequence of I|1U|) is that Aui can at 
best be adjusted in steps of \5ui s — Su>i\. With such steps, 



the frequency mismatch Aw can always be made to be- 
long to the interval [— \ s "*~ Sui \ ; l^-fo^l j canno ^ b e 
brought closer to zero. This minimal Alo corresponds to 
the mode pair chosen by the OPO. By considering the 
worst case, the maximal frequency mismatch that can be 
reached is 



Acij max — 



\Su) s - 8u)i\ 



(11) 



which agrees with the general expression of the detun- 
ing given in |l6j |. As the OPO is pulled away from an 
exact cavity resonance by increasing the cavity length, 
the frequency mismatch will increase up to the maximal 
value given by 111(1 . at which point there will be an op- 
erating mode with a lower mismatch to which the OPO 
will switch. This phenomenon is well known as mode 
hopping [THE!- 

Before 1)110 can be compared with the bound found for 
the Hopf bifurcation in Sec. [H] it has to be expressed 
in the same units. The detunings used in equations 
are normalized so that half-height width of the cavity 
resonance is 2. Since the cavity finesse T is defined so 
that the half- height width is Slo/J 7 , we have 



A, 
A, 



2T t 
2Ti 



Alo s 

' Suj s 
Au t 
SuJi 



(12a) 
(12b) 



As mentioned before, A s and Aj are not independent 
in the stationary regime but obey the simple relation 
A s = Aj as a consequence of energy and photon number 
conservation Taking into account that Alo = Alo s + 
Awj, one obtains from i|12|) : 



A., = A,; = 



2T,T,Al 



T 1 5lo s + F s 8lo 1 



(13) 



This expression was derived assuming lo s ~ u>i , so that we 
should fix T s ~ Ti for consistency. The maximal value 
of the frequency mismatch authorized by mode selection, 
as given by equation l|ll|) . can thus be rewritten as 



A M - T 

"max — •'t 



\Su s - 5ujj\ 
Slo s + 5u)i 



(14) 



Before we can discuss whether the Hopf bifurcation can 
be observed in typical experiments, we have to reformu- 
late 1 )140 in terms of the experimental configuration. 

Given that 



2ttc 



2{L C av + {l7>s,i l)y 



(15) 



where L cav is the geometric length of the cavity, l c is 
the length of the nonlinear crystal, n s ,i the indices of the 
signal (resp. idler) fields and c is the celerity of light in 
the vacuum, 1)140 can be rewritten as 



\Sn\l c 
[L] 



(16) 
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where 



[L] = 2(L cav + (?±±* - l)l c 



is the average optical path for the signal and the idler 
fields over one round trip in the cavity and 



Sn 



In the limit case of a monolithic OPO (L cav = l c ), 
expression <|16|) leads to the remarkably simple expression 



A M IT = 

"max/'' s 



2\Sr 



(17) 



showing the key role played by cristal birefringence. Ex- 
pression (|17|) yields an absolute upper bound for i|16|) 
since cavity length obviously cannot be smaller than crys- 
tal length. 

Note that since Sn <C n s ,i for a standard crystal, ex- 
pressions l|16[) and H17|l ensure that the small-detuning 
hypothesis of the mean-field approximation is fulfilled for 
the signal and idler fields. 

We are now in a position to obtain a simple criterion 
to determine whether the Hopf bifurcation can be ob- 
served in a given configuration at infinite pump power. 
Obviously, we must have 



.m 



A • < A 

mm ■ "max 



(18) 



which, using (jSJ and (|16|) . and expressing the pump cav- 
ity decay rate as 



7 : 



^1 

T ' 

P 



translates into 




(19) 



which is the main result of our paper. Again, l|19|l sim- 
plifies in the monolithic case. A noteworthy feature of 
inequality l|19fl is the asymmetry in the dependences with 
respect to the pump and signal finesses. It is easily seen 
that increasing pump finesse is much more effective to 
have l(T§|l satisfied. Assuming T p = 50 and T s = 500, 
doubling the pump finesse decreases the left hand side 
of (|19|) by 46% while doubling the signal finesse only 
does so by 4%. 



B. Numerical estimates 

To get a better understanding of whether the crite- 
rion obtained in Sec IIII Al is easily satisfied or not, we 
now compute numerical estimates for typical experimen- 
tal configurations. In previous experiments [l3 , ll-'l IbH 
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FIG. 2: Solid lines : values of maximal detuning authorized 
by mode selection when cavity length is varied between 1.5 
cm and 10 cm. Dashed lines : minimal detuning permitted by 
the conditions of the Hopf bifurcation for two different values 
of 7 : 5 and 10. 



0,|2(j, we used a KTP crystal of length l c = 15 mm cut 
for type-II phase-matching, with an extraordinary index 
n e — 1.75 and an ordinary index n a — 1.83. The crys- 
tal is enclosed in a Fabry-Perot cavity delimited by two 
spherical mirrors with a radius of curvature of 5 cm. The 
cavity finesses for the signal and pump fields are around 
500 and 50, respectively. 

The most easily adjustable parameter is the geomet- 
rical cavity length, which is bounded from below by the 
crystal length (1.5 cm) and from above by the concentric 
condition (^10 cm). Figure [21 shows the evolution of the 
maximal detuning A^f ax with cavity length. The main 
feature is that A^f ax decreases monotonically with cavity 
length, as is easily seen in expression <|16|1 . Thus the most 
favorable situation will always be obtained in the "mono- 
lithic" configuration where cavity length equals crystal 
length. In practice, this configuration cannot be reached 
when spherical mirrors are used because of the size of the 
crystal mount, but the shortest feasible cavity length (3 
cm in our experiments) should be sought. 

In our experiment, 7 = T s /T p — 10, from which we 
can compute the value A^ in = \J 7(7 + 2) of the minimal 
detuning necessary to obtain the Hopf bifurcation. This 
value does not depend on cavity length and is represented 
in Fig.Elby the upper dashed horizontal line. We see that 
there is hardly a configuration where A^ in < A^ ax . In 
this configuration, the Hopf bifurcation cannot be ob- 
served regardless of pump power because mode hopping 
prevents to reach sufficiently high detunings. 

In order to lower A^ in , assume now that pump finesse 
is increased to 100 so that 7 decreases to 5. This cor- 
responds to the lower horizontal dashed line in Fig. [21 
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There is now a range of cavity lengths where A min < 
A^f ax as shown by the black area between the two curves. 

If T s and T p are increased simultaneously keeping 7 
constant, A^f av will increase while A^- remains un- 

7 nidx mm 

changed, and the instability region will also widen. 



C. Conclusion 

In this section we have taken into account the fact 
that the signal frequency detuning is not a fixed parame- 
ter but is determined by the operating frequency chosen 
by the OPO so as to minimize this detuning. As a re- 
sult the signal frequency detuning is bounded from above 
by the mode hopping phenomenon. Since signal detun- 
ings at which periodic behaviors appear are bounded 
from below, the Hopf bifurcation can only be observed if 
A^ in < A^f ax . Our analysis has shown that even in the 
ideal case where infinite pump power is available, mode 
hopping can prevent from reaching detunings sufficiently 
large to observe the Hopf bifurcation of the monomode 
model. 

The numerical estimates obtained for the configuration 
used in our previous experiments explain why we could 
not observe the Hopf bifurcation in this setup. In order 
to obtain instabilities, one should take into account that 
A^ in depends on the ratio of signal and pump finesses 
while A* J ax is proportional to signal finesse. More gen- 
erally, the finesses should satisfy inequality (|19|l . As will 
be discussed in Sec. [V] the theoretically most favorable 
configurations are extremely difficult to build experimen- 
tally. In all cases, cavity should be made as short as 
possible. 

It remains to be checked that in cases where the insta- 
bility can be observed, it persists when pump power is 
limited. We do so in Section llVl 



IV. NUMERICAL INVESTIGATIONS AT 
FINITE PUMP POWER 

So far, our analysis has assumed infinite pump power, 
which has allowed us to obtain simple analytical formu- 
las such as <|19[) to decide whether the Hopf bifurcation 
is screened by mode hopping or not. We now have to 
determine how good this approximation is in the real-life 
situation where only finite pump power is available and 
if conclusions drawn from our theoretical analysis remain 
relevant. 

At finite pump power, whether the Hopf bifurcation 
can be observed no longer depends only on inequal- 
ity (0} but also on whether the Hopf threshold © can 
be reached given available pump power. Because expres- 
sion JSJ is much more complicated than inequality , we 
restrict ourselves here to a numerical exploration of the 
detuning ranges where periodic behavior is found and 
compare our results with predictions from the infinite 
pump power analysis. This exploration is carried out 



for several values of the cavity finesses, other parame- 
ters b eing chosen so as to match our previous experi- 
ments 1 1 la . 

EE Q EE As we shall see, it will allow 
us to conclude that pump power is not a limiting factor 
with commonly available pump lasers and that the cri- 
teria derived in the infinite pump power analysis remain 
relevant. 

A few general observations are in order before we 
present our numerical results. The relevant criterion is 
whether the minimum detuning A^ in at which Hopf bi- 
furcation occurs is smaller than the larger detuning A^ ax 
allowed by mode hopping. Thus it is interesting to com- 
ment on how these bounds evolve when pump power is 
limited. 

Regarding onset of periodic behavior, it should be re- 
called that inequality Q holds regardless of pump power, 
that equality can be achieved only for infinite pump 
power and that otherwise the Hopf threshold (JSJ) yields a 
more stringent condition than Q on the detunings. As 
a result, the instability regions are systematically shifted 
towards higher values of the detunings, leading to an in- 
crease of A^ in compared to expression © . 

As for the maximal value A^ ax of the detuning al- 
lowed by mode hopping [Expr. I|14[l]. it does not depend 
on pump power as it is obtained by considering the fre- 
quency combs of the cavity resonances. However a limi- 
tation that has to be taken into account at finite pump 
power is that if the OPO is below parametric emission 
threshold at A s = A^xj then the latter bound certainly 
cannot be achieved and the actual bound will be lower 
(the stationary ON state must exist for the bifurcation 
to occur). Using the expression for parametric thresh- 
old Pjl. the expression for the maximal signal detuning 
value becomes 



AiZ = min A£, 



E 2 



1 



A 2 

p 



- 1 



(20) 



where E is the pump parameter and A p is the pump 
detuning. The compatibility between the values of the 
minimal detuning at which periodic behavior occurs and 
of the maximal detuning at which stationary paramet- 
ric emission occurs is now more difficult to analyze be- 
cause both depend on pump detuning. However, it is 
easy to see that as pump power is decreased, the former 
can only increase and the latter only decrease so that the 
Hopf bifurcation is necessarily harder to observe in the 
finite pump power case. We now assess by how much by 
carrying out numerical simulations. 

The phase diagrams in Fig.|21shows in black the regions 
in the (A p ,A s ) parameter plane where periodic behav- 
ior is found for different values of the cavity finesses, the 
values of the remaining parameters being fixed so as to 
match our experiments |12l Ha. 1 1 ll Ulx |20| . For each pa- 
rameter set, the dynamical regime is classified as periodic 
when 

• Inequality is satisfied, and 

• Eh(A p , A s ) < E max {T v ,Ts) 



7 



where E ma ^(J- p , J- s ) is the maximum pump parameter 
corresponding to the pump power available in our exper- 
imental conditions (4 W) [23. Since the unstable zones 
are enclosed inside the regions where parametric emission 
occur and the instability regions at infinite pump power, 
the boundaries of these regions, which can be computed 
analytically, are also shown in Fig. [31 so that we can es- 
timate how well they approximate the numerical results. 
In order to make meaningful comparisons between setups 
corresponding to different values of the finesses, the do- 
mains of variations of the pump and signal detunings are 
chosen so that A S:P E [— J- atP /(10ir), J- s , p /(10n)]. This 
corresponds in each case to the same variation of the 
physical cavity length (AL = A/407r) and thus there is 
no difference in the scans from an experimental point 
of view. Since |A SjP | <c 2J-" SjP , this also ensures that 
the small-detuning hypothesis of the mean-field model is 
satisfied. 

Fig. EJa) shows that for parameter values correspond- 
ing to our experimental setup, periodic behavior does not 
occur, explaining why we do not observe the Hopf bifur- 
cation in this setup. It also shows that the boundaries 
of the unstable region at infinite pump power are located 
outside the central band where A s e [— A* ax , A*f ax ] and 
thus that mode hopping prevents the bifurcation from be- 
ing observed in this configuration even if infinite pump 
power was available. The proximity of the two curves 
corresponding to parametric threshold and bifurcation at 
infinite pump power probably explains why no instability 
can be observed. 

In contrast with Fig. |3{a), Figs. |3Jb) to Old) display 
larger and larger instability zones as cavity finesses are 
increased. While in Fig. E2b) the intersection of the un- 
stable region with the central band of allowed signal de- 
tunings is very small, it becomes sizable in Fig. |2{d). 
More precisely, the fraction of the central band occupied 
by unstable regions is 1.6% in Fig.^b), 3.5% in Fig.[3tc) 
and 17% in Fig. |3{d). These numbers are meaningful as 
probability estimates if we assume that there is no cor- 
relation between the signal and pump detunings on av- 
erage, i.e., that the entire allowed band may be explored 
over several experiments. It should be noted that during 
a scan of cavity length through a single resonance of the 
pump, pump and signal detuning will vary jointly as illus- 
trated in Fig.|3Ja) and that this may affect the probabil- 
ity of hitting the unstable zone. However, it is expected 
that there is no relation between paths followed in the 
(A p , A s ) plane for different pump resonances so that no 
part of the allowed band should remain inaccessible. 

Examination of Figs. Et a )-(d) shows that although in- 
creasing finesses globally makes it easier to observe the 
Hopf bifurcation, finesses for the pump and signal fields 
have different influences. The twofold increase in T v be- 
tween Figs. Ola) and (b) clearly modifies the phase dia- 
gram much more than the twofold increase in T s between 
Figs. OJb) and (c). Similarly, increasing T v from 100 
to 500 is critical to have a significant probability of ob- 
serving the Hopf bifurcation. This is consistent with the 
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FIG. 3: Phase diagrams in the (A P ,A S ) plane showing 
the instability regions as black areas for different values of 
the finesses of the pump {Tp) and of the signal (J- s ) field: 
(a) T v = 50, T s = 500, (b) T v = 100, T s = 500, (c) 
T p = 100, T s = 1000, (d) T p = 500, T s = 1000. Solid curves 
delimit areas above parametric emission threshold. Dashed 
lines indicate where |gj is an equality and delimit unstable 
areas in the infinite pump power approximation. Dotted hor- 
izontal lines indicate the maximal detuning value allowed by 
mode hops. In (a), the dotted and dashed line shows a pos- 
sible path followed by detunings as cavity length is scanned. 
The discontinuity occurring as the path reaches the horizontal 
dotted line corresponds to a mode hop. 
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discussion of criterion i|19|) in Sec. IIIII which showed that 
increasing pump finesse was much more effective than in- 
creasing signal finesse. However increasing pump finesse 
is extremely difficult from a practical point of view, as 
we discuss in Sec.[V] Interestingly, it was similarly noted 
in Ref. [ToT ] that increasing pump finesse while keeping 
signal finesse allowed one to observe an hexagonal trans- 
verse pattern at lower signal detunings. 

An important conclusion that can be drawn from the 
numerical exploration summarized in Figs. Ola)-(d) is 
that except in Fig. [3Ja), the theoretical analysis at in- 
finite pump power provides a very good approximation 
of the finite pump power case, all the better as finesses 
are higher and as observing the Hopf bifurcation becomes 
more plausible. Indeed the instability zones are tightly 
delimited by the boundaries obtained at infinite pump 
power and by the parametric threshold line except for 
larger pump detunings. In particular, the agreement is 
excellent at the tip of the unstable zones, near the point 
of minimal signal detuning. This makes us confident that 
for pump lasers currently available, the criterion obtained 
in 119|) is effective in assessing the probability of occur- 
rence of the Hopf bifurcation. It also shows that pump 
power is not a limiting factor as when there is a nonempty 
intersection between the unstable zone and the central 
band of allowed detunings, its area only marginally in- 
creases from the finite to the infinite pump case, mostly 
in regions far from the minimum detuning value. 

To conclude, numerical calculations at finite pump 
power show that while the Hopf bifurcation becomes 
harder to observe than in the infinite pump power case, 
the mode hopping phenomenon remains the main limit- 
ing factor by limiting the range of values signal detun- 
ing can take. In order to obtain experimental evidence 
of the bifurcation, building a dedicated setup would be 
much more effective than increasing pump power. How- 
ever, there are some experimental difficulties in doing so, 
which we discuss in Section Ivl 



V. EXPERIMENTAL CONSIDERATIONS 

The theoretical analysis at infinite pump power of 
Sec. IIIII and the numerical computations at finite pump 
power of Sec. IIVI have suggested that the Hopf bifurca- 
tions of the monomode mean-field model might be ob- 
servable in some configurations. In this Section, we dis- 
cuss the feasibility of an experimental setup specially de- 
signed for evidencing the bifurcation. 

Some parameters are easily optimized using the find- 
ings of Sec. IIIII and IIVI As discussed in Sec. IIII Bl cav- 
ity should be as short as possible and in this respect a 
monolithic configuration would be optimal. Then, crys- 
tal birefringence should be chosen as large as possible 
[see Eq. (|17|l ]. However this may not be an option as 
brifrigence is primarily used to achieve phase-matching 
in the desired operating range. Of course, pump power 
should be as large as possible, but our analysis has shown 



that there were modest returns in increasing it much 
beyond that offered by currently available pump lasers. 
The two parameters left for setup optimization are then 
the signal and pump cavity finesses. Criterion l|19D indi- 
cates that those finesses, and especially the pump finesse, 
should be taken as large as possible. 

However, cavities of very high finesse are critical to 
align and are very sensitive to fluctuations. For example, 
a change of 5L of the cavity length will induce a variation 
of 5 A = AT8L/X of the detuning, where A is the optical 
wavelength. For A = 1064 nm and J- = 1000, a fluctua- 
tion SL = 5 A of the cavity length will induce a variation 
of the detuning SA = 2, corresponding to full-width of 
the resonance. The length of the cavity must thus be ad- 
justed carefully and maintained constant to within less 
than a few angstroms. 

This stability problem is all the more critical when 
thermal effects are taken into account. In our configura- 
tion, it has been shown that thermal effects can induced 
complex dynamical behaviors 0,0] during which cavity 
length can be spontaneously swept by a few nanometers. 
Such instabilities were observed even at incident pump 
powers around 500 mW 12) , which for a pump finesse of 
45 corresponds to an intracavity pump power of about 
6.5 W. For 4 W of incident pump power and a finesse 
of 500, intracavity pump power would reach ~640 W 
at resonance and thermal effects would then be a hun- 
dred times as strong. For illustration purposes, let us 
recall that the only chaotic regime reported so far in an 
TROPO was observed in a situation where cavity length 
could not be made stationary 0] , although the configu- 
ration was standard. 

Another adverse influence of pump absorption in the 
crystal is that it puts a limit on the highest finesse achiev- 
able. In our setup, for example, pump absorbtion in the 
1.5cm-long crystal is 2% cm -1 which implies that pump 
finesse cannot be higher than about 100 even if perfectly 
reflecting mirrors were used. Except for crystals with ex- 
ceptional low absorption, high values of the pump cavity 
finesses such as 500 appear to be completely unrealistic. 

Apart from the previous remarks which hold for all 
systems featuring an absorbing material enclosed in a 
resonant cavity, specific restrictions on the finesses are to 
be considered for a TROPO, inside which several fields 
must be simultaneously resonant. Generally cavity fi- 
nesses in an experimental TROPO setup are chosen so 
that T v <C T s . It ensures that numerous resonances of 
the signal field are found inside the wider resonance of 
the pump field, in other terms that several mode hops 
occur before OPO falls below threshold. This allows one 
to only consider coincidences of signal and idler cavity 
modes. For very high cavity finesses of the pump, the 
mode with smallest signal detuning might very well be 
below threshold because of a large pump detuning. In 
fact, the analysis of part IIIII should then be reworked 
and reformulated in terms of coincidences of three reso- 
nances instead of two. It is not certain that this more 
complicated problem admits solutions. 
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To conclude, an OPO designed in order to observe the 
Hopf bifurcation should have a short cavity (e.g., mono- 
lithic setup), high birefringence and pumping rate. The 
pump cavity finesse should be as high as possible. Under 
these conditions, however, an OPO would be very diffi- 
cult to operate and to stabilize on an operating point. 
Even then, the numerical analysis of Sec. IfVI indicates 
that unstable zones would remain small compared to the 
range of operating parameters. This casts serious doubts 
on the experimental observability of the Hopf bifurcation 
of the monomodc mean-field model. 



VI. CONCLUSION 

In this paper we have raised the issue that contrary to 
what is commonly assumed in theoretical studies of the 
dynamics of optical parametric oscillators, signal detun- 
ing is not a fixed parameter but a dynamical variable 
whose value results from a complex selection process. 
Thus, detuning limitation due to mode hopping may af- 
fect the occurrence of dynamical instabilities occurring 
at high values of the detuning. As a simple example, 
we analyzed under which conditions the Hopf bifurcation 
leading to periodic behavior in the monomode mean-field 
model of a triply resonant OPO can be observed. 

We first showed that signal detuning must reach a min- 
imal value for the bifurcation to occur, even at infinite 
pump power. We then described how the mechanism of 
mode hopping restraints the values signal detuning can 
take and we gave an analytical expression of the maximal 
detuning allowed by mode hopping. By comparing the 
two bounds so obtained, we showed that they are incom- 
patible in many configurations, showing that finite pump 
power is not necessarily the limiting factor for observing 
this instability. This was confirmed in numerical compu- 
tations of instability domains carried out at a finite pump 
power corresponding to our experimental configuration, 
in which we observed that most of the unstable domains 
are inaccessible because of mode hopping. We found that 
pump finesse is a critical parameter that should be made 
as large as possible in order to observe the Hopf bifur- 
cation. However experimental setups designed for this 
purpose should most certainly be extremely difficult to 
operate. 



For the sake of simplicity, our derivation of the max- 
imal detuning allowed by mode hopping relied on two 
main hypotheses. One was that signal and idler have 
frequencies close to 1:1 ratio but experience different in- 
dices, as in our type-II experiments. Generalization to 
other frequency ratios in the type-II configuration should 
be straightforward and is not expected to modify the 
conclusions of the present work. Type-I OPOs near fre- 
quency degeneracy have more complex tuning properties 
and should be analyzed separately. The other hypothe- 
sis was that parametric gain can be considered constant 
along a frequency domain containing many longitudinal 
modes so that the lowest possible detuning is always se- 
lected. If that is not the case, then values of signal de- 
tuning higher than predicted may be obtained and our 
analysis would have to be adapted. However, it does not 
seem likely that our conclusions would be modified a lot. 

An open question is whether conclusions drawn for 
the Hopf bifurcation of the monomode mean-field model 
also hold for other theoretical predictions of temporal 
or pattern-forming instabilities in doubly or triply opti- 
cal parametric oscillators. Indeed, such instabilities have 
usually been predicted for relatively high values of the 
signal detuning and/or of the pumping rate (see, e.g., 
Refs. HH H H, 0, Il3)- Thus, their observability may 
very well be also affected by detuning limitation due to 
mode hopping rather than by finitcness of the pump 
power. A detailed comparison of instability thresholds 
in detuning space with the maximum values allowed by 
mode hopping is in order for these theoretical predictions 
so that their relevance for experimental systems can be 
assessed. To conclude, mode hops are not only a nuisance 
for stabilizing and tuning OPOs but also for using them 
as tools to study complex dynamical behavior. 
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